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The electronic bandstructure and the Fermi surfaces of ferromagnetic CeRh 3 B 2 are calculated by using 
FLAPW and LSDA+ U method. As assuming several kinds of the ground state to describe the Af electronic 
state, we propose a fully orbital- and spin-polarized state |/- = 0, s x = 1 /2> as the ground state, instead 
of the conventional LS -coupled CEF ground state, generally expected in typical 4/ compounds. This is 
supported by the fact that both the observed magnetic moment and the observed dHvA frequencies are 
C*) , well explained by the calculated electronic structure and the Fermi surfaces. The unconventional ground 

D ■ state is stabilized by the strong 4/-4/ direct mixing between the neighbored Ce atoms along the extremely 

P-h ' small distance along the c-axis in the hexagonal crystal cell. 
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1. Introduction 

The trivalent cerium compounds show the versatile fea- 
tures such as the heavy fermion behavior, 1 ' 2 ' the anisotropic 
superconductivity, 1 ' the complicated magnetic phase, 3 ' 4 ' and 
the multipole order. 5 ' These phenomena demonstrate both the 
localized and itinerant behavior of Af electrons. Among the 
ferromagnetic cerium compounds with other nonmagnetic el- 
ements, ternary cerium boride CeRfr^ has the highest Curie 
temperature Tq =120 K, 6 ' while the usual magnetically 
ordering temperature is only 1-10 K in other cerium com- 
pounds. In contrast to the highest 7*c, it is notable that the 
saturated magnetic moment in CeRh3B2 is remarkably small. 
The magnetic measurement done by Galatanu ef.a/. 7 ' at 2K 
has shown that the saturated magnetization shows the strong 
anisotropy in the hexagonal crystal structure: 0.45 IfiB along 
[1010] ; 0.447yu B along [1 120] ; 0.04 jig along [0001] direction 
so that they have concluded that the easy axis is the [1010] di- 
rection. This small values of the magnetic moment (^0.45^) 
are considerably smaller than the value (»1.0 yUg/Ce) 8 ' of gen- 
eral ferromagnetic Ce compounds. In order to explain the un- 
usual ferromagnetic behavior, a Ce 4/ itinerant ferromagnetic 
model, 9 ' a Rh Ad itinerant ferromagnetic model 10 ' and Ce 4/ 
localized ferromagnetic models 7, n ~ 16 ' have been proposed so 
far. 

In the Af localized model, the ground state of the Ce 4/ 
electron is described by the crystalline electric field (CEF) 
splitting and the relatively stronger spin-orbit splitting. In the 
hexagonal symmetry, the six-fold degenerate \j = 5/2) level 
splits into three doublets, where the ground states is described 
by \j z = ±1/2). However, this scheme fails in the explanation 
of the observed small magnetic moment because the doublet 
ground state \j z = +1/2) is expected to cause very large satu- 
rated moment of 1.6 hb in the basal plane. Therefore, the Af 
localized model with considering a CEF ground state hardly 
explain the observed small magnetic moment. 
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On the other hand, the de Haas-van Alphen (dHvA) 
measurement, observed for LaRli3B2 and ferromagnetic 
CeRli3B2, 17 ' strongly supports the 4/ localized model. This 
observation has suggested that the contribution of the Af elec- 
trons to the Fermi surfaces is small (i.e. the Af electron is 
well localized) because the topology of the Fermi surfaces 
in CeRli3B2 is similar to that of LaRli3B2, which is consis- 
tent with the theoretical result 18 ' calculated by a full potential 
LAPW (FLAPW) method. Therefore, the localized Af model 
can be applied to explain the Fermi surfaces, however, as men- 
tioned above, there is inconsistency in the localized model in 
the magnetic moment. 

In this paper, we report the theoretical study of the elec- 
tronic structure of CeRli3B2, aiming at clarifying the property 
of the localized Af electron. First we will explain the detail 
of the method of the calculations, then we will discuss the 
ground state of the Af electron, which leads to the particular 
magnetism and the Fermi surfaces. 

2. Computational Details 

We have performed a bandstructure calculation within the 
local density approximation (LDA) in the density-functional 
framework using the full potential linearized augmented 
plane-wave (FLAPW) formalism. 19 ' We used TSPACE and 
KANSAI-99 program codes for this calculation. 

The scalar relativistic effects are taken into account for all 
electrons and the spin-orbit interactions are included for the 
valence electrons inside the Muffin-tin spheres as in a sec- 
ond variational procedure. 20 ' In the perturbation Hamiltonian 
for the spin-orbit coupling (SOC), the radial potential of the 
spin-offdiagonal element is substituted as an average of the 
potential for the spin-up and down state. 21 ' The SOC Hamil- 
tonian term is unitary transformed by using the spin-rotation 
matrix, which defines the quantization axis so that the cal- 
culated electronic state is dependent on the direction of the 
magnetic moment. The net orbital moment is induced by the 
SOC which breaks the time-reversal symmetry. 

In order to treat carefully the strong correlated Af state, 
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the effective L(S)DA+f/ potential 22 ' is applied to improve 
the LSDA calculation. For the Af electron system, where the 
spin-orbit coupling is involved, the density matrix is described 
as a spin- and orbital-dependent (i.e. 14 x 14) matrix, con- 
sidering the spin-flip element. 23 ' The effective potential is in- 
troduced into the calculation as a second variation, together 
with the SOC Hamiltonian. In this procedure, although the 
spin- and orbital-dependent density matrix converges self- 
consistently with the wavefunctions, the resulted Af state is 
strongly dependent on the initial density matrix with a certain 
value of U. In other words, the ground state of the 4/-state is 
determined a posteriori. Based on the same method, a meta- 
magnetic transition of Fermi surfaces of CeRu2Si2 has been 
successfully investigated. 24 ' In the following sections, we will 
first report the bare LSDA result and later the LSDA+t/ re- 
sults with considering the several Af ground states. 
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Fig. 1. Crystal structure of CeRfi3B2; (a) The hexagonal structure contain- 
ing three unit cells and (b) the convention unit cell of the base-centered or- 
thorhombic structure with the lower symmetry in the ferromagnetic state. 
The Cartesian coordinates used in the calculation are also shown. 



CeRli3B2 crystallizes in a hexagonal CeCo3B2 crystal 
structure (Z)' space group), shown in Fig. 1(a). It should be 
emphasized that the hexagonal unit cell is used only for the 
calculation of the non-magnetic state whereas the less sym- 
metrical orthorhombic unit cell is used for the ferromagnetic 
state. This is because the symmetry of the electronic state is 
broken by the SOC, which couples the freedom of spin and 
orbital, so that the reduced symmetry is dependent on the 
direction of the collinear magnetic moment. As taking into 
account of the previously reported magnetic measurement, 71 
which shows that the magnetic moment is oriented in the basal 



plane with the small in-plane magnetic anisotropy, we have 
adopted the base-centered orthorhombic unit cell (Z)^? space 
group, shown in Fig. 1(b)), where the magnetic moment is as- 
sumed to be aligned along x ([1210] or equivalently [1 120] in 
the hexagonal structure) axis, and also another setting where 
the magnetic moment is aligned along the y ([1010]) axis is 
considered, for the comparison. 

The experimentally measured structural parameters 25 ' are 
used in the calculation : a = 5.474 A and c = 3.085 A (Note 
that b — a in the hexagonal cell and b — V2a in the or- 
thorhombic cell). Muffin-tin radii are set as 0.2761a, 0.2192a 
and 0.1761a for Ce, Rh and B sites respectively. The core 
electrons (Xe-core except 5s 2 and 5p 6 for Ce, Kr-core ex- 
cept 4p 6 for Rh, He-core for B) are calculated inside the MT 
spheres in each self-consistent step. The LAPW basis func- 
tions are truncated at |k + G,| = 5.81 x 27r/a, corresponding 
to 409 LAPW functions at the T point. The sampling 150 k- 
points (divided by 6, 6 and 10) are uniformly distributed in 
the irreducible l/8th of the orthorhombic Brillouin zone. 

3. LSDA result 

Before we proceed to the LSDA+J7 calculation, let us focus 
on the bare LSDA result, where the Af state is treated as an 
itinerant state. Considering the ferromagnetic state, the spin 
magnetization is set to be parallel to the x axis ([1 120] direc- 
tion in the hexagonal structure as in Fig 1(a)). In this study, 
we use the definition that the spin and orbital quantization is 
described with respect to z axis ([0001] axis in the hexagonal 
structure). Therefore, the magnetic quantum number m of the 
spherical harmonics of the LAPW basis corresponds to the 
z-projected orbital moment (L). 

Figure 2 shows the calculated bandstructure for CeRli3B2 
nonmagnetic and ferromagnetic states within the L(S)DA 
scheme. They show the similar aspects except for the tiny 
spin-splitting of the / bands near the Fermi energy in the 
ferromagnetic bandstructure, where the up- and down- spin 
states are mixed by SOC. In the energy region which is shown 
in the figure, the bandstructure mainly consists of Rh-d states 
below the Fermi energy and the Ce-o" states above the Fermi 
energy. Figure 3 shows the density of states for the nonmag- 
netic state. The main components of density of states at the 
Fermi energy are Ce-/ (48%), Rh-c/ (18%) and Ce-d (17%) 
states. Due to the small polarization of both the spin and the 
orbital states, the value of the net magnetic moment along x 
axis is negligibly small as 0.07 fig in the ferromagnetic state. 
Therefore, the bare LSDA calculation, based on the itinerant 
/ model, fails to reproduce the observed magnetic moment. 
In addition, the calculated Fermi surfaces (not shown) are 
completely different from what expected from the observed 
dHvA observation, 17 ' which is similar to those of prototype 
LaRli3B2. This is because the calculated / bands have itiner- 
ant nature, as crossing the Fermi level. 

In order to improve the LSDA result, we introduce the ef- 
fective LSDA+{/ potential into the Ce-4/ state. In this pro- 
cedure, the effective potential shifts down the occupied /' 
level by — U/2 and shifts up the unoccupied / 13 levels by 
+U/2 with respect to the original LSDA levels. The occu- 
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Fig. 2. Bandstructure of CeRh3B2 for (a) the nonmagnetic state and (b) the 
ferromagnetic state calculated in the bare LSDA scheme. The first Bril- 
louin zone for (c) the nonmagnetic hexagonal lattice and (d) the ferromag- 
netic orthorhombic lattice is also shown. In (b) and (d), £' (2/3, 0, 0) and 
A' (2/3, 0, 1/2) points and YP (2/3, 0, k z ) line are newly defined. 
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Fig. 3. An LDA result of the density of states for nonmagnetic CeRh3B2. 
Partial density of states of Ce-/ (highlighted), Ce-d (dashed line) are Rh-rf 
(dotted line) states are indicated. 



pied f 1 level in the LSDA bandstructure is labeled as "fo" 
in Fig. 2. This state has been already discussed by Takega- 
hara et al. within the APW calculations. 26 ' They have pointed 
out that due to the strong ddcr and ffcr mixing along the ex- 
tremely short Ce chain along the z axis, the bottom of the 
Ce-4/o (= F30) and 5do (= F20) bands are shifted downward 
largely with respect to other / and d bands and show the large 
dispersion in energy. This mechanism is also explained by 
the CEF point charge model: the fo wave function tends to 
extend to the nearest neighbor positively charged Ce 3+ ions. 
In our LSDA result, the spin-degenerated fo state remains in 
the ground state even in the ferromagnetic phase: the fo state 
strongly hybridizes with the conduction state and does not 
show the spin-split. In order to obtain the appropriate ferro- 
magnetic state, we must remove the degeneracy of fo spin 
states. 

4. LSDA+U result 

4.1 Choice of the 4/ ground state 

According to the above discussion, we have chosen the 
fully spin- and orbital- polarized state f^ p = \l~ = 0, s x = 1 /2) 
as the 4/ ground state in the following LSDA+tV calcula- 
tion. This choice may conflict with the conventional idea, 
where the Af ground state is described with the strongly spin- 
orbit split \j) state with the comparatively smaller CEF split- 
ting. In CeRli3B2, however, our choice is reinforced by the 
assumption that the fo state is stabilized by the larger CEF 
splitting than the SOC splitting, due to the extremely short c 
length. The \j = 5/2) state and \j = 7/2) state hybridize well 
and compose the |/) and \s) states separately, in spite of the 
\j) states. Therefore we treated the \l z ) and \s x ) states sepa- 
rately to construct the density matrix. The spin state is set 
as fully polarized along x (or y) axis according to the direc- 
tion of the observed magnetic moment. After constructing the 
initial density matrix where one electron occupies f " p state, 
the effective potential shifts down the f^ 1 ' band and shift up 
the other / bands. Then the density matrix is calculated self- 
consistently so that the final ground state is determined. 

4.2 Bandstructure and the magnetic moment 

Figures 4 and 5 show the calculated bandstructure and the 
density of states, respectively, in the LSDA+t/ scheme where 
U is set as 0.3 Ry. The Ce-d and / levels are artificially shifted 
upward by 0.1 3Ry and 0.1 2Ry respectively. This modifica- 
tion follows the previous FLAPW calculation for LaRh3B2, 18) 
where the La-d and / levels are shifted by 0.1 Ry and 0.2Ry 
to obtain the proper Fermi surfaces. Such modification is nec- 
essary to obtain the proper Fermi surfaces, as reported in 
LaB 6 , 27) YbAl 3 , 28) YNi 2 B 2 C, 29 ' 30) LuNi 2 B 2 C 31) study, espe- 
cially when a boron atom is involved in the system. 

The LSDA+t/ effective potential changes drastically the 
bandstructure with respect to the LSDA result. Near below 
the Fermi energy, the f^ p band is fully occupied as show- 
ing the large dispersion, which reflects the symmetry of T30 
wavefunction and the strong /-/ mixing. On the other hand, 
the spin-split d^ and di ovm bands are fully unoccupied above 
the Fermi level. Therefore only Kh-d bands cross the Fermi 
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Fig. 4. Calculated bandstructure of the ferromagnetic CeRfi3B2 in the 
LSDA+[/ scheme (a) with the [lOlO] magnetization and (b) with the 
[1120] magnetization. The parameter U is set as ().3Ry. 
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Fig. 5. Calculated density of states for ferromagnetic CeRh3B2 in the 
LSDA+t/ scheme with the [1010] magnetization. The notation is the same 
as in Fig. 3. 



level, showing the large spin-splitting at the Z point and along 
the YP axis. Comparing the bandstructure with the magneti- 
zation along x (Fig. 4 (a)) and along y (Fig. 4 (b)) axis, the 
dependence of the bandstructure on the direction of magneti- 
zation is considerably small. The finally obtained f 1 ground 
state has not only the main component of |/- = 0, s x =1/2) 
state (as the initial set) but also the small components of 
\l z = ±1, s x = 1/2) state because the SOC term has the non- 
zero matrix element between the different \L) states by ±1. In 
consequence, the obtained ground state is described as 

{a \l z = 0) + b\l z = ±1)} \s x = 1/2) 

where a = 0.98 and b = -0.15 are obtained. The mixing be- 
tween different \L) state produces the orbital magnetic mo- 



ment along the x axis. As listed in Table I, with the com- 
parison to the experimental result, the calculated spin mo- 
ment (s x = 0.98/^b; s y = 0.90fiB) of Ce-/ state is almost 
fully-polarized and the orbital moment (l x = -0.85/ig; l y = 
-0.81/ig) cancels the spin moment whereas Ce-d and Rh-ii 
states don't have any significant magnetic moment. Therefore, 
the net magnetic moment is only 0. 1 1 (0. 16)iub/Cc along x (y) 
axis. This reduced moment is the key which explains the ob- 
served small saturated magnetic moment of 0.45 ps- 

4.3 Spin-split Fermi surfaces 
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Fig. 6. Fermi surfaces of the ferromagnetic CeRh.3B2 : (a) hole surface 
from the 60th band, (b) hole surface from the 61st band, (c) hole surface 
from the 62nd band, (d) hole surface from the 63rd band, (e) electron sur- 
face from the 64th band, (f) electron surface from the 65th band, (g) elec- 
tron surface from the 66th band, (h) electron surface from the 67th band 
and (i) electron surface from the 68th band with the [1 120] magnetization: 
(J) electron surface from the 67th band band with the [1010] magnetization. 



In contrast to the nonmagnetic LaRh^Bi, which has the five 
Fermi surfaces, 18 ' we have obtained the spin-split nine Fermi 
surfaces in ferromagnetic CeRli3B2 as shown in Fig. 6. The 
spin-splitting of the Fermi surfaces, which mainly consists of 
Rh-ii state, is due to the magnetic contribution from the spin- 
and orbital-polarized / ground state. We have also checked 
here the dependence of the Fermi surfaces on the magnetiza- 
tion axis. As changing the magnetization axis, it is found that 
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Table I. The orbital and spin magnetic moment (jib) of ferromagnetic CeRti3B2: The first three lines show the experimental result. The fourth and fifth lines 
show the calculated results with the magnetization is aligned along [1120] and [1010] axis, respectively. 
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Fig. 7. Angular dependence of the calculated dHvA frequencies of the fer- 
romagnetic CeRfi3B2 with (a) the [1120] magnetization and (b) the [lOlO] 
magnetization. The direction of the field is denoted with the orthorhombic 
lattice vector, together with the corresponding direction in the hexagonal 
lattice at the bottom of the figure. 



Table II. Experimentally observed dHvA frequencies F exp and calculated 
dHvA frequencies F ca i c in units of 10 3 T, cyclotron masses of m exp & m ca i c 
in unit of the free electron mass and mass enhancement factor A. The mag- 
netic field is applied along the [0001] and [lOlO] direction. The experi- 
mental value is taken from Ref. 17 ' 
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only the 67th electron surface (shown in Fig. 6 (h) and (j)) 
give the slightly different shape. In fact, when the magnetic 
field is applied to the ferromagnetic system and the magnetic 
moment is flipped, SOC changes the shape of Fermi surfaces 
according to the direction of the magnetic moment. 

The angular dependence of the dHvA frequency derived 
from the Fermi surfaces is shown in Fig. 7. Although the 
shape of the Fermi surfaces can change progressively accord- 
ing to the applied magnetic field, we used rigid Fermi surfaces 
independent on it. Recalling that the magnetic moment is ori- 
ented in the basal plane, it may be sufficient to consider the 
Fermi surfaces only with the [1120] and [1010] magnetiza- 
tion. Between the two cases, the significant difference is due 
to the 67th electron surface whereas the rest show the similar 
frequencies and angle dependence. The dHvA branches a, flu 
($2, 82 and s are named by the experimental results whereas the 
branches a, b, c and d are named here. The dHvA frequencies 
and the electron mass are compared with the experimental re- 
sults, 17 ' as listed in Table II. In the rest of this subsections, we 
closely look the comparison and discuss the spin-split effect 
on the Fermi surfaces. 

The 60th and 61st ellipsoidal hole surfaces correspond to 
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the 30th surface of LaRh3B2 where the splitting is relatively 
small. The two dHvA branches named a show the satisfactory 
agreement with the experimental result. 17 ' 

The 62nd and 63rd hole surfaces correspond to the 3 1 st sur- 
face of LaRhsBa. The closed spherical Fermi surfaces named 
P\ shows similar aspect of the spin-splitting to the case of a 
surfaces, however, only one of the branches has been detected 
in the measurement. On the other hand, the large spin-splitting 
between the 62nd surface and the 63rd surface is essential to 
reproduce the experimental Fermi surfaces. The splitting sep- 
arates the closed 62nd surfaces (02) and the connected 63rd 
surface (62) forming a ring along k z axis. It results in the fact 
that the dHvA branch is observed when the field is applied 
along [0001], whereas the 62 branch is not observed. This re- 
sult is at variance with the experimental expectation by Okubo 
et.al., ll) where they have attributed the 62 branch to the 64th- 
67th sheet-shape electron surfaces. 

The 64th-67th electron surfaces have typical quasi-one- 
dimensional flat shapes which correspond to the 32nd and 
33rd surfaces of LaRh3B2. The 67th surfaces has the closed 
orbits, a, b and c, however, the corresponding dHvA branches 
show the different angular dependence from the observed 62 
branch, therefore we insist again on that 62 is the spin-split 
counterpart of the branch. 

The 68th ellipsoidal Fermi surface centered at the Z(0, 0, 
1/2) point correspond to the 34th surface of LaRh3B2. The 
counterpart of the surface disappears due to the large spin 
splitting. The s branch show the similar angular dependence 
with the experimental result, however the value of the cross 
section is underestimated here. 

To summarize the above discussion, we have obtained the 
proper Fermi surfaces, which explain well the experimentally 
observed dHvA branches. The angular dependence of a, 
P2, 62, s branches is substantially in good agreement with the 
experiment, while the cross section of £2 and s branches are 
underestimated in the calculation (cfr. Table. II). 

The missing branches in the calculation are j\ and 72 
branches, observed with H || [1120] and with H || [0001], 
respectively. In the experimental study, Okubo et.al. has at- 
tributed both of them to the 64th-67th flat surfaces. As com- 
paring the angular dependence, the calculated b orbital may 
correspond to the observed j\ branch whereas we cannot re- 
produce the 72 branch. Although the angular dependence of 
the observed 72 branch implies that there is a large closed or- 
bit perpendicular to k z axis, it is hard to construct such a large 
orbit from these flat surfaces without changing the volume of 
the electron surfaces, but the volumes for electrons and holes 
must be compensated. Therefore we guess the 72 branch orig- 
inates from another surface, such as 63rd hole surface. The 
hypothesis may be valid if one assumes larger spin-splitting 
between the 62nd and 63rd surfaces so as to enhance the vol- 
ume inside the 63rd hole surface. In such situation, the cross 
sectional area of hole orbit 62 is increased and, at the same 
time, the area of an hole orbit inside the ring-shaped sur- 
face is decreased. The hole orbit with the small area causes 
the dHvA branch, which property may be consistent with the 
experimentally observed 72 branch. In this context, the cross 



sectional area of d orbital is also enhanced and therefore the d 
branch may be shifted upward, out of the observable range. It 
explains why the calculated d branch has not been observed. 

More serious problem is that the spin-splitting of fi\ 
branches have not been experimentally observed (i.e. only one 
of the pair is observed). Since f}\ surfaces have the similar 
band character (mainly Ce-d and B-p character) and the cross 
sectional area to those of a surfaces, the spin-splitting effect 
in dHvA frequency is supposed to be in proportion to the cy- 
clotron mass, which is observed approximately twice as large 
as that of a orbital. Therefore it is supposed that fi\ branches 
must show the spin-splitting twice as large as a branches, as 
indeed shown in our calculation result. So far, there is no clear 
reason to explain why only the spin-splitting of a branch has 
been observed but not the one of f}\ branch. 

As we end this subsection, it should be emphasized again 
that the spin-splitting is important to discuss the observed 
dHvA result for ferromagnetic CeRli3B2. The occupied f? p 
band doesn't cross the Fermi level so that the topology of the 
Fermi surfaces is similar to that of LaRlv^, however, the hy- 
bridization from Jq P band, which is very close to Fermi level, 
causes the spin-splitting for the conduction bands. In order to 
discuss the effect, let us have a careful look at the bandstruc- 
ture. Figure 8 shows the bandstructure of the ferromagnetic 



CeRh.B, 



o 10-49% 
Ce/component q 5f>|0o% 



i.- • 




Fig. 8. Bandstructure of the ferromagnetic CeRfi3B2. The Ce-f component 
is shown as round symbols in the bands. See the text for the details. 



CeRli3B2 with the [1120] magnetization within the enlarged 
energy scale. The bands which have large Ce-f component is 
indicated with round symbols, which size stands for the ratio 
of the component. It is shown that the cf interaction is large 
around the Z and T points and along YP axis, then the con- 
ductive Rh-d band is largely spin-split. The splitting in the 
vicinity of the Z point corresponds to the s surface splitting 
and the splitting along YP axis corresponds to the ^2-62 sur- 
face splitting. 

4.4 The effect of the choice of ground state 

In the above LSDA+f/ calculation, we have chosen the 
/q state as the ground state, instead of the conventional CEF 
state. In order to check the validity of the consideration, an- 
other LSDA+jy calculation was performed, where the ground 
state is chosen as a hexagonal CEF state \j = 5/2, j : = 1/2} 
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state. In this choice, it is assumed that the Kramers pair 
\j — 5/2, j. — ±1/2) of the CEF level in hexagonal symme- 
try is spin-split so that one of the pair is fully occupied. The 
spin magnetization is set as parallel to z axis ([0001 ] direction 
in the hexagonal symmetry). 




r ' ' r * a s hs'l r am " lk p ii 

Fig. 9. Calculated bandstructure of the ferromagnetic CeRfi3B2 in the 
LSDA+{/ scheme with the [0001] magnetization. The parameter U is set 
as 0.3Ry. 

Figure 9 shows the calculated bandstructure of the ferro- 
magnetic CeRh3B2 with \j = 5/2, j z = 1/2) ground state with 
U = 0.3 Ry. In this calculation, the Ce-d and / levels are 
artificially shifted upward by 0.13 Ry and 0.08 Ry respec- 
tively as trying to fit the Fermi surfaces to the experimental 
dHvA result. The calculated orbital and spin magnetic mo- 
ment is 0.45 and -0.18 /-Ib along [0001] axis, which is com- 
pared with the expectation value of the CEF level as 0.57 
and -0.14 fj.%, respectively. The topology of Fermi surfaces 
(not shown) reflects the hexagonal symmetry and the dHvA 
frequency shows different angular dependence from the prior 
LSDA+C/ result; the and 62 branches are widely split, but 
the spin-splitting at a branch is not seen. This is fatal to 
reproduce the experimental Fermi surfaces so that the prior 
LSDA+J7 calculation with f^ p ground state is favored. Be- 
sides, the magnetic moment along [0001] axis of this ground 
state is totally different from the experimental result. 

5. Conclusion 

Here we briefly review the characteristic of Fermi surface 
at magnetic 4/ electron system as focusing on the localized 
vs itinerant behavior. It is known that spin-split Fermi sur- 
faces are observed at magnetic Pr and Nd compounds, as ex- 
plained by localized / model. For example, filled skutterudite 
NdFe4Pi2 shows clear spin-pairs of dHvA branches which 
shows small spin-splitting, as compared to LaFe4Pi2. 34 ' This 
shows the localized character of 4/-electrons in Nd and the 
existence of small but sizable c-f interaction. By contrast, 
Ce compounds, due to the rather stronger c-f interaction, are 
classified into several groups as following: 

1) A localized model is valid at CeAl 2 and CeB 6 , 35 ~ 37 ' 
where the Fermi surfaces are quite similar to non-/ reference 
system, LaAF and LaBf,, as indicating that the 4f electron is 
essentially localized and contribute little to the formation of 
the Fermi surfaces. 



1') However, CeSb • ' at ferromagnetic phase shows the 
quite large spin-splitting of dHv A branches (e.g. fiu fh, /?3 and 
/J4), where one cannot easily assign these branches in pairs 
due to the various angular dependence, because anisotropic 
hybridization based on the local /-electron affects the Fermi 
surface topology. 

2) At CeSn3, the 4f electrons are considered as itinerant 
electrons; The Fermi surfaces are well explained by band cal- 
culation, where the / electrons are treated in the same manner 
as usual s, p, d conduction electrons. 40 ' 

Recently, heavy fermion superconductivity without spacial 
inversion symmetry has been observed at CePtiSi and related 
compounds, CeTXi (r=Rh, Ir, Co; Z=Si, Ge). 41) It is theoret- 
ically expected that each Fermi surface of magnetic CePtiSi, 
a, /3 and y is spin-split by the Rashba-like antisymmetric spin- 
orbit coupling. 42 ' Besides, even nonmagnetic reference sys- 
tem, LarGe3 (r=Fe, Co, Rh, Ir) shows the split Fermi sur- 
faces. 43 ' Then no further spin-split is observed in this system, 
even if it is magnetically ordered. Such the large spin-orbit 
band splitting due to the lack of the inversion symmetry is 
considered as the origin of the particular symmetry of the su- 
perconducting gap node. 

CeRu2Si2 shows metamagnetic transition from the non- 
magnetic ground state to the magnetic state at H m ~7.7T, 
where both the occupied Ce-/ 1 state and the magnetic mo- 
ment concomitantly change. 44 ' A recent LSDA+tV calcula- 
tion of the Fermi surfaces with and without externally applied 
magnetic field explains well the observed dHvA branches at 
both states, showing that some of dHvA branches are largely 
spin-split in the applied magnetic field. 24 ' The transition is at- 
tributed to the change of Fermi surfaces from 2) to 1'). 

At Ce7Tn5 (r=Co, Rh, Ir) series, known as heavy fermion 
superconductors, FLAPW band calculations with itinerant 
4/ model well explain the Fermi surfaces of Celrlns and 
CeCoIn5, 45 ' 46 ' whereas the Fermi surfaces of CeRhlns shows 
the high deviation from them, which implies the localized na- 
ture of the 4/ electrons. 47 ' Therefore, the /electron behavior 
changes a lot even in the similar compounds of CeTIns series. 

Here we note that CeRli3B2 has similar aspect to CeSb 
case. The above theoretical results compared with the ex- 
perimental indicate that only LSDA+t/ calculation with /"'' 
ground state can explain the experimentally observed Fermi 
surfaces at ferromagnetic CeRb.3B2. It is reasonable that due 
to the non-appropriate treatment of electron correlation, a bare 
LSDA calculation fails to reproduce the observed Fermi sur- 
faces as well as the expected magnetic moment. We empha- 
size that at this system, the spin-splitting of the conduction 
band is strongly affected by the c-f interaction so that it is 
not obvious to find the spin-split pairs of experimentally ob- 
served dHvA branches. Therefore LSDA+t/ calculations are 
necessary to assign the dHvA branches to each Fermi sur- 
face, where we showed in this paper the ground state of Ce- 
/ state and the direction of the magnetic moment changes 
the shape of Fermi surfaces. Along this line, we can unam- 
biguously specify the electronic state of 4/ electron system 
using by LSDA+J7 method and comparing with the dHv A re- 
sult. The validity of the f^ p ground state can be confirmed 
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by other means of experiment: Indeed, a preliminary result of 
polarization-dependent photoelectron spectra recently shows 
the evident peak of f° state below the Fermi energy. 48 ' 

In summary of this paper, the Fermi surfaces of ferromag- 
netic CeRh3B2 is theoretically investigated. We showed that 
the spin-splitting at the conduction band is caused by the mag- 
netic contribution from fully spin-polarized Ce-/o band. By 
considering the /q ground state and the effect to the conduc- 
tion bands, we explained well the observed small magnetic 
moment as well as the experimental dHvA result. The large 
band dispersion of fi state due to the /-/ direct mixing be- 
tween Ce sites may give a hint on the origin of the High Tq at 
CeRh 3 B 2 . 
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